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of equations is decoupled in the sense that the correlations of In this paper we address nonper- 
turbative aspects of the analytic theory of hydrodynamic turbulence. Of paramount importance 
for this theory are the "fusion rules" that describe the asymptotic properties of n-point correlation 
functions when some of the coordinates tend toward one other. We first derive here, on the basis 
of the Navier- Stokes equations, a set of fusion rules for correlations of velocity differences when all 
the separation are in the inertial interval. Using this set of fusion rules we consider the standard 
hierarchy of equations relating the n-th order correlations (originating from the viscous term in the 
Navier-Stokes equations) to n + I'th order (originating from the nonlinear term) and demonstrate 
that for fully unfused correlations the viscous term is negligible. Consequently the hierarchic chain 
of equations is decoupled in the sense that the correlations of n + I'th order satisfy a homogeneous 
equation that may exhibit anomalous scaling solutions. Using the same hierarchy of equations when 
some separations go to zero we derive a second set of fusion rules for correlations with differences 
in the viscous range. The latter includes gradient fields. We demonstrate that every n'th order 
correlation function of velocity differences Tn{Ri, R2, ■ ■ ■) exhibits its own cross-over length r;„ to 
dissipative behavior as a function of, say, R\. This length depends on n and on the remaining 
separations R2, R3, ■ ■ ■■ When all these separations are of the same order R this length scales like 
VniR) ~ rt{R/LY'^ with x„ = (C„ — C,„+i + Cs - C2)/(2 — C2), with Cn being the scaling exponent of 
the n'th order structure function. We derive a class of exact scaling relations bridging the exponents 
of correlations of gradient fields to the exponents (,„ of the n'th order structure functions. One of 
these relations is the well known "bridge relation" for the scaling exponent of dissipation fluctuations 



I. INTRODUCTION 

In a recent series of papers |,|,|,| we developed an 
approach to the theory of the universal scahng proper- 
ties that characterize the statistical invariants of fully de- 
veloped turbulent flows. The main technical tool of this 
approach was renormalized perturbation theory which al- 
lows, after exact resummations, to obtain insights about 
the non-perturbative properties of the universal statistics 
of turbulence. The main nonperturbative results which 
were obtained are (i) Hydrodynamic interactions are lo- 
cal in physical-space and in scale-space once the sweeping 
effects are removed, (ii) The inner, viscous scale which 
appears perturbatively as an Ultra- Violet cutoff leads 
non-perturbatively to the anomalous scaling of correla- 
tions involving velocity gradients, (iii) The outer scale 
of turbulence appears in the theory as a renormalization 
scale after infinite resummations of a renormalized per- 
turbation series that converges order by order. The aim 
of this paper is to present new nonperturbative results 
that together with the previous findings begin to assem- 
ble to a solid structure of the theory of hydrodynamic 



turbulence. Brief presentations of some of the new re- 
sults were offered in |||,§,0,D . 

The issue under study is not new; decades of ex- 
perimental and theoretical attention (see for example 
1^,^ 12 l^,|l^) have been devoted to two types of simul- 
taneous correlation functions; the first type includes the 
structure functions of velocity differences, 
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(1.2) 



where (...) stands for a suitably defined ensemble aver- 
age. A second type of correlations include gradients of 
the velocity field. An important example is the rate e(r, t) 
at which energy is dissipated into heat due to viscous 
damping. This rate is roughly v\'Vu{r,t)\'^. An often- 
studied simultaneous correlation function of e{r,t) — 
e{r, t) — e is 



K,,{R) = (i{r + R)iir)) 



(1.3) 



Within the famous Kolmogorov 1941 (K41) approach 
|]l5| to turbulence one predicts that for R in the in- 
ertial range, i.e rj ^ R ^ L, Sn{R) is {eRY^^^ and 
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K^^{R) ~ v'^t^^^R~^/^ . Here L and rj are respectively 
the outer scale of turbulence and the Kolmogorov vis- 
cous cutoff, e is the mean energy flux per unit time per 
unit mass. 

Experimental measurements show that in some aspects 
Kolmogorov was remarkably close to the truth. The ma- 
jor aspect of his predictions, that the statistical quantities 
depend on the length scale R as power laws, is corrobo- 
rated by experiments. On the other hand, the predicted 
exponents seem not to be exactly realized. For example, 
the correlation K^g{R) decays according to a power law, 



K^,{R) oc R- 



(1.4) 



for rj R <^ L, with fi having a numerical value of 
0.2 — 0.3 instead of 8/3 jl^. The structure functions also 
behave as power laws. 



Sn{R) OC i?^" , 



(1.5) 



but the numerical values of Cn deviate progressively from 
n/3 when n increases |p^ , p^ . Something fundamental 
seems to be missing. 

The first positive answer within our approach appeared 
in the context of correlation functions of gradient fields 
rather than of the velocity differences [p|jlq]. As indi- 
cated above the discrepancy in the exponent fi between 
the experimental value and the naive expectation is large. 
This gave hope that the explanation must lie close to the 
surface. Indeed, considering the perturbative scheme for 
K^^{R) leads immediately to the discovery of logarith- 
mic ultraviolet divergences with some ultraviolet viscous 
scale fj acting as the renormalization scale. The summa- 
tion of this infinite series results in a factor [R/fjY^ with 
some anomalous exponent A which is, generally speak- 
ing, of the order of unity. Additional calculations lead 
to the exact result A = 2 — C2 and to the understanding 
that this result means that the renormalized perturba- 
tion series for the structure functions diverges in the limit 
L — > cxD like {L/RY". The anomalous exponents 5„ are 
the deviations of the exponents of 5„(i?) from their K41 
value. This is a very delicate and important point. It is 
not simple to see why a series each of whose terms con- 
verges when L — > 00 still diverges as a whole. One can 
understand this as in re-summing the series of a desired 
function one finds an inhomogeneous equation whose in- 
homogeneous solutions are indeed L-independent. How- 
ever, the equation possesses also homogeneous solutions 
which are inherently nonperturbative in nature and may 
have anomalous scaling. Homogeneous solutions must be 
matched with the boundary conditions, and this is the 
way that the outer scale appears in the theory. An im- 
portant remaining step in the theoretical development is 
to understand how to compute the anomalous exponents 

The divergence of the entire perturbation series for 
Sn{R) with L — > 00 forces us to seek a nonperturbative 
handle on the theory. One fact that may potentially lead 
to a nonperturbative control is the existence of a global 



balance between energy input and dissipation. This may 
be turned into a nonperturbative constraint on each n- 
th order structure function Using the Navier-Stokes 
equations one derives the set of equations of motion 



dSniR,t) 
dt 



D^{R,t) = J„iR,t) 



(1.6) 



where Z)„ and J„ stem from the nonlinear and the vis- 
cous terms in the Navier Stokes equations respectively. 
In the stationary state the time derivative vanishes and 
one has the balance equation 



Dn{R) = Jn{R) 



(1.7) 



The evaluation of Dn{R) does not pose conceptual dif- 
ficulties. It was shown that it is of the order of 
dSn+i/dR. The evaluation of Jn{R) poses some deli- 
cate issues, but it was shown that it is proportional to 
Sn{R)/S2{R)- This result is of course the famous hierar- 
chic chain of equations of which many closure attempts 
have failed to provide the desired anomalous scaling so- 
lutions. 

At the heart of this paper we off er a way out of this 
conundrum. The problem with Eq.(l/7) is that it refers 
to a "fully fused" quantity, and this fact leads to difficul- 
ties. What is meant by this is that the primitive object 
of the statistical theory of turbulence is not the structure 
function Sn{R), but rather the "fully unfused" n-rank 
tensor correlation function of velocity differences: 



F„(ri,r'i;r2,r' 



1) ' 2, ' 2i ■ ■ ■ i '^nj '^n) 

(w{ri,r[)w{r2,r'^) . . .w{rn,r'J) . (1.8) 



In this quantity all the coordinates are distinct. The 
structure function Sn (R) is obtained by fusing all the co- 
ordinates Ti into one point r, and all the coordinates r'^ 
into another point r + R. In this process of fusion one 
crosses the dissipative scale. To control this process we 
need to formulate the so-called "fusion rules" for turbu- 
lence, which govern the properties of n-point correlation 
functions when two or more coordinates "fuse" together. 
One cannot fuse blindly. 

The main fundamental results of this paper are as fol- 
lows: 

1. Fusion rules in the inertial interval. The explana- 
tion of how fusion works was briefly given in 
We begin this paper in Sections II and III with a 
detailed derivation of the fusion rules, considering 
all the possible relative geometries of the sets of 
fused and unfused points. 

2. The homogeneous equation for the {n + 1)- order 
correlation functions: the fusion rules allow us to 
consider in Section IV and V the generalized bal- 
ance equation (4.7) for the fully unfused correlation 
functions. We demonstrate that for these quan- 
tities the hierarchic chain of equations decouples, 
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since now the dissipative term vanishes when the 
viscosity v goes to zero. In the case of Eq.(|l.7|) the 
term J„ contains a correlation function with many 
fused coordinates. As a resuh (and cf. Section VI) 
it remains finite in the hmit ^ 0. In the present 
case of fully unfused quantities we get in the limit 
i/ — > a homogeneous equation involving only cor- 
relations of order n + 1. The full analysis of this 
equation is not easy, and it deserves a considerable 
amount of attention in the future. 

3. The existence of "viscous scaling functions". Hav- 
ing obtained the generalized balance equation we 
can now study precisely when, as a function of one 
separation distance (say \r[—ri\), the fully unfused 
correlation function crosses over to dissipative be- 
havior. It turns out that the scale at which this 
happens depends on the rank n and on the remain- 
ing inertial range distances. The dependence on the 
inertial range distances is characterized by anoma- 
lous exponents, and we compute them exactly in 
terms of the anomalous exponents C„. This is done 
in Section VI. 

4. Exact bridge relations. With a precise evaluation of 
the viscous scaling functions and the fusion rules in 
hand we can consider correlation functions of gradi- 
ent fields, and relate their scaling exponents to the 
anomalous exponents Cn- The result of this exer- 
cise, which is presented in Section VII, is an infinite 
set of bridge relations. We consider correlations of 
the type 

^p"^ = {e{xi)e{x2) . . . e{xp)w{ri,r\). . .w{r,,, r'^)) 

(1.9) 



where i? is a typical separation between any pair, 
and we are interested in the scaling relations be- 
tween the exponents fiii and the exponents C„. 

(2) 

Note that /Xg in this notation is the well studied 
[0,0 exponent of dissipation fluctuation which is 
denoted ji. On the basis of the Navier-Stokes equa- 
tions and the fusion rules we establish rigorously 
that 



(p) 



n+3p 



(1.10) 



In particular we offer a solid derivation of the 
phenomenologically guessed bridge relation [p|Jl4| 
A* = 2-C6. 

In Section VIII we summarize the paper and indicate the 
direction of the road ahead. 



II. FUSION RULES FOR THE INERTIAL 
INTERVAL 

In this section we discuss the fusion rules that arise 
when some of the coordinates of an n-point correlation 



function tend toward one other. We begin by introducing 
the set of correlation functions that are required by the 
analysis. 



A. Correlation functions and related quantities 



The fully unfused correlation function Fn was intro- 
duced in Eq.(l.S). Other statistical quantities of interest 
have the same number of velocity differences but they 
depend on a smaller number of coordinates. In the lan- 
guage developed below we will refer to such quantities as 
partially fused correlation functions. For example if all 
the velocity differences are evaluated with respect to the 
same point we define an n-rank tensor 



S„(ro|ri,r2. . .r„) = (u;(ro, ri)u;(ro, r2) . . . i(;(ro, r„)) . 

(2.1) 

Note that in an inhomogeneous ensemble this quantity 
depends explicitly on rg. Obviously, Fn depends on 
many more coordinates than Sn- It is useful sometimes, 
therefore, to represent correlation functions Fn in terms 
of the partly fused quantities Sn. This can always be 
done: represent w{rj,r'j) as w{rj,ro) + w{ro,r'j), and 
by substitution find the needed relation. It can also be 
useful to choose Tq as the center of mass of all the coor- 
dinates in Fn. We will return to such relations when we 
need them. 

One can fuse more points together. For example we 
can have only one separation vector distance, and form 
the n-rank tensor 



SniR) = Sn{ro\r) = Sniro\r,r . . .r) , 



(2.2) 



where i? = r — rg. Equat ion s ( |2.lD and (2.2) can be un- 
derstood as the limit of ( |l.8[ ) when pairs o f co ordi nate s 
fuse together. For example, in going from (1.8) to (2.1) 
we fuse all the unprimed coordinates, and keep all the 
primed coordinates different. To get ( |2.2D we also fuse 
all the primed coordinates. 

Note that one can also consider another type of fusion, 
in which pairs of coordinates rj and r'j which are as- 
sociated with the velocity difference w{rj,r'j) coalesce. 
We will refer to the fusion of coordinates with a velocity 
difference as a fusion of type A, and a fusion of coordi- 
nates that are not associated with a velocity difference as 
a fusion of type B. The properties of these two limiting 
procedures are different, and they have to be studied sep- 
arately. It is commonly assumed that the limit of type B 
exists, and the limit remains finite when ^ 0. This is 
equivalent to the statement that S'„(rg|ri, r2 . . . r„) and 
SniR) are independent of Re for asymptotically large 
Re. Although this property looks innocent, it is a deep 
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statement that in general needs to be supported by the- 
ory. In the language of diagrammatic perturbation the- 
ory this property is a consequence of the locality of every 
diagram in the perturbation expansion for these quanti- 
ties; this property was discussed in detail in In ||] 
it was explained that there exists a mechanism for the 
appearance of the viscous scale in the statistical theory. 
However, this mechanism operates only when one corre- 
lates gradients of the velocity fields rather than the ve- 
locity differences themselves. Thus the above statement 
is equivalent to the assumption that there are no addi- 
tional non-perturbative mechanisms for the appearance 
of the viscous scales in the theory of correlations of veloc- 
ity differences. Additional objects that appear naturally 
in the theory are correlations whose tensorial nature is of 
lower rank, including scalars and vectors. The first one is 
the scalar quantity which is appropriate for even orders 
of iS'„. To keep in mind its scalar nature we denoted it 

in Q as S2m (R) and define it as 

S2m{R)^{\w{ro,r)\^"') , R = r-ro. (2.3) 

o 

The quantity 5*2171 (R) is analytic. For odd order struc- 
ture functions we introduce a vector object S2m+iiR) 
according to 



^2m+l{R) 



{w^{ro,r)\w{ro,r)\'"') 



(2.4) 



Here and below we will use Greek indices to indicate vec- 
tor and tensor components, and Roman indices to indi- 
cate the order of the quantity. The placement of indices 
as subscripts or superscripts has no special meaning, and 
is chosen for convenience. 

In isotropic turbulence the vector S2m+i{R) can only 
be oriented along R. This allows us the introduction of 
a scalar quantity S2m+iiR) which depends on the mag- 
nitude of R: 



R 

S2m+l{R) = -^S2m+l{R) 



(2.5) 



In later sections we use the objects defined here to study 
the fusion rules and their consequences. 



B. Basic Properties of Correlation Functions in 
Systems with Flux Equilibrium 

The systems that we discuss are driven on a character- 
istic scale referred to as the outer scale L. This driving 
can be achieved by either a time dependent low frequency 
"stirring force" or by specifying given values of m at a set 
of "boundary" points with a characteristic separation L 
away from our observation points rg, Tq, ri, etc. The sys- 
tem suffers dissipation due to viscosity, and in the inviscid 
limit the kinetic energy is conserved. In three dimensions 
we deal with a "direct" energy cascade in which the in- 
take of energy on the scale L is balanced by dissipation 
on a small scale rj ^ L. 



We invoke two fundamental assumptions [|5| : 

1. Scale invariance: all the correlation functions are ho- 
mogeneous functions of their arguments in the core of the 
inertial interval 77 ^ jr^ — ro| ^ £: 

F„ ( Ar 1 , Ar 1 ; . . . ; Ar„ , Ar ^ ) = A'^" F„ (r 1 , r'^ ; . . . ; r„ , r ^ ) , 

(2.6) 

where C„ are scaling exponents. It is obvious that all the 
other types of correlation and structure functions that 
emanate from this most general quantity and which were 
detailed in the previous subsection share the same scal- 
ing exponents. Accordingly ^„ is also the usual scaling 
exponent of the n-th order structure function (2.3). 

2. Universality of the scaling exponents: this means that 
we can fix an arbitrary set of velocity differences on the 
scale of L, and the scaling exponents of the correlation 
functions will not depend on the precise choice of the 
L-scale motions. Mathematically this is expressed as the 
following property of the conditional average: if for i < n 
\ri — r^l ^ L, and for i > n, \r'^ — ri\ ^ L, then 

{w{ri,r[)w(r2,r'2) . . . w{rn,r'J\w{rn+i,r'^+i) 
■ ■ ■ w{r„+N, <+7v)) = Fniri,r[; . . . ; r„, r'J 
X ^nMirn+i,r',^+i; ■ ■ .;r,,+N,r'„+iy) (2.7) 

The precise meaning of the universality assumption is 
that the functions Fn have the same scaling exponents 
as Fn in the inertial interval. They may be different 
functions in the inertial interval and in particular they 
may differ in their crossover to viscous behavior. Their 
(different) crossover scales may depend on the large scale 
motions that were fixed in the conditional average. 



C. Fusion rules for the fusion of p pairs of points in 
velocity differences (type A) 

1. General case: 2 < p < n — 2 

The first set of fusion rules that we derive concerns 
Fn when p pairs of points ri,r[ . . . rp,rp, {p < n) which 
involve p velocity differences tend to some point rg, see 
Fig. II C 1 . Here we exclude the special cases p — 1 and 

p = : 



— 1 in which the leading scaling contribution van- 
ishes by symmetry. We will consider these special cases 
in the next subsections. 

Consider the situation in which the typical separation 
of all the coordinates ri, . . . ,rp from Tq is r, whereas 
all the other separations remain much larger, say of the 
order oi r <^ R <^ L. We will show, on the basis of 
assumptions 1 and 2, that 

F„(ri,r;;...;r„,0 (2.8) 
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where Fp is a tensor of rank p associated with the 
first p tensor indices of The (n — p)-rank tensor 

*n,p('"p+ij ■'"p+i; ■ • • ; '"nj rJJ is a homogeneous function 
with a scahng exponent ^„ — Cpi ^ind is associated with 
the remaini ng n — p indices of Fn ■ The derivation of the 



fusion rule (2^) follows from Baycs' theorem. We write 



F,,n,r;...r„,<)=/<i»,r,„,<„).. 

X {w{ri,r[),w{r2,r'^) . ..w{rp,r'p) ■w(rp+i, r^+i) 

X w{rp+2,r'p+2) ■ ■ ■ w{rn,r'J) , (2.9) 

where V[w{rp^i,r'p^i) . . . w{rn,r'^)] is the probability 
to see the tensor w{rp+i,r'p_^i) . . . u>(r„, rj^). Next note 
the consequence of assumption 2: the scaling laws of the 
correlation functions at scale r are the same independent 
of whether we force the system on the scale L ^ r or 
on the scale R ^ r. The conditional average in (|2.9|) is 



proportional to Fp, and hence (2.8) 



FIG. II C 1 The geometry of fusion of type A. Lines connect- 
ing points indicate velocity differences across that distance. In 
this example there are three velocity differences across small 
separations (of the order of r) and two velocity differences 
across large separations, of the order of R. 



The first impression is that (2.8) means statistical inde- 



pendence of the small scale motion from the large scales. 
This is not so. The difference is that the assumption of 
statistical independence would lead to Fn = FpFn-p- 
Indeed, the first factor Fp has the same order of mag- 



nitude and the same exponent as Fp in (2.8). However 



the factor F^-p has a scaling exponent Cn-p rather than 
Cn ~Cp which is the exponent of ^n,p- When all the large 
separation are of the same order i?, *n,p Sn{R) / Sp{R) 
which in the case of multi-scaling is much larger than 
Fn-p ^ Sn-piR)- We thus understand that the fusion 
rules in fact demonstrate the existence of a very special 
statistical dependence of the small scales on the large 
scales. This dependence stems physically from the exis- 
tence of a direct energy flux from large to small scales. 
We will see that it can lead to a totally unconventional 
scaling structure of the theory. We should stress that 
these fusion rules were derived from first principles for 
p = 2 in Navier-Stokes turbulence 0| and for passive 
scalar advection for any p fia] . 



2. Fusion of one pair of points in a velocity difference 

As we mentioned the case p = 1 in which the velocity 
difference across a small scale r = \ri — r[\ is correlated 
with n — 1 velocity differences across larger distances re- 
quires a special attention. A naive application of the 



fusion rule (2.8) results in Fn oc Fi{ri,r[) which van- 
ishes due to space homogeneity. In order to evaluate the 
leading non- vanishing term we expand u(r[) in a Taylor 
series around u{ri), 

u{r[) = u(ri) + ViM(ri) • {r[ - n) + . . . . (2.10) 

Using this we can write 

lim F„(ri,ri; . . . ;r„,r^) 

= (r'l - ri) • Vi (w(ri)w;(r2, r^) . . . w{rn, <)) . 

The correlation function in this formula depends an all 
the separation distances, and the gradient with respect 
to ri picks up contributions from all the differences 
rj — ri. Therefore the gradient can be evaluated as the 
inverse of the smallest of these, |Vi| ^ l/Rmin where 
i?min = mini{|ri — ri|}. For r <^ i?min this leads to the 
evaluation 



F„(ri,ri; . . . ;r„,r^) 



Rn 



-SniR) 



(2.11) 



This formula creates an immediate worry about the situ- 
ation in which iJmin becomes of the order of r or smaller. 
We will analyze this situation in Subsection A. Now 
we will consider the next special case of fusion, when 
p — n — 1. 



3. Fusion rule for one large separation distance 

Consider F„(ri, r'^; r2, r'j; . . . ; r„, r^) with all the co- 
ordinates being nearby except for ri which is far away, a 
distance R from the remaining 2n — 1 coordinates which 
are all within a ball of radius r. Since we assume that the 
flow is isotropic, the tensor Fn is a (generally reducible) 
representation of the rotation group. Moreover, in the 
situation discussed here the isotropization of the small 
scales (of the order of r) with respect to the direction of 
R = r'i~ ri leads to a direct product structure 

i^,r""-""(r-i,r'i;r2,r;,;...;r„,r:j 

= C:^Fr ■""(^,r-2,r;,;...;r„,r;) . (2.12) 

The constraint of incompressibility, which is always true, 
can be written as 



^1 n 



'(ri,r;;r2,r^;...;r„,0 = . (2.13) 



Applying the divergence operator to ( 2.12| ) we find the 
following form of the incompressibility constraint. 



C\2 + R^ ] F„"^ ■""(i?,r2,r^;...;r„,,r:. 







(2.14) 



This equation has two solutions: (i) C = 0, (ii) Fn oc 
1/R^. However according to the general rule F„ must be 
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proportional to R^"~'^"-'^. Thus the first solution is real- 
ized, the contribution which is oc vanishes due 
to the incompressibility constraint. Therefore we need 
to consider the next non-vanishing evaluation which is 
i?-indepcndent : 



;r„,r;)«r^". (2.15) 



This result is interesting. It means that in the case con- 
sidered here in which all the separations of order r except 
one (which is of order R, r] ^ r R L) the scaling ex- 
ponent of the correlation function Fn is fully determined 
by small r-scale fluctuations. 

In addition to these cases one needs some special ge- 
ometries of fusions of type A in which the general evalua- 
tion is inapplicable. These cases are treated in Appendix 
A and referred to as needed. 



D. Rules for Fusions of Type B 



FIG. [ID The geometry of fusion of type B. Lines connect- 
ing points indicate velocity differences across that distance. In 
this example there are three velocity differences across large 
separations, but three coordinates within a ball of small ra- 
dius r. 



1. Coalescence of two points (p = 2) 

Fusions of type B refer to situations in which there 
are p coordinates within a ball of small radius r, but no 
two coor dinates that belong to a velocity difference, see 
The simplest case is p 



Fig II D. The simplest case is p = 2 in which \ri — 
is much smaller than any other separation. Consider the 
velocity differences w{ri,r[) and w{r2,r'2). They can 
be reexpressed as 



w{ri,r[) = w{ri,ra) +w{ro,r[) 
w{r2,r'2) = w{r^,r'2) +w{r2,ro) 



(2.16) 



where tq = [ri + r2]/2. This allows us to write 

{w{r,,r[)w{r2y2){wr-^) = {w{r[,r,)w{r'2,ro){wr-^) 
-(u;(ri,ro)u;(r2,ro){w}""^) - {w{r[,ro)w{r2,r„){w}''''^) 
+ (-u;(ri,ro)-u;(r2,ro)M"-2) . (2.17) 

Here we used the short-hand notation {w}"^^ to denote 
the remaining product of n — 2 velocity differences. 

The first term on the RHS is independent of the small 
separation, and is a homogeneous function of the large 
separations with a scaling exponent The next two 
terms contain one velocity difference across small sep- 
arations which according to the discussion of Fig. is 
proportional to r^^ . The last term has two velocity differ- 
ences across a small separation, and is also proportional 



to the same factor. This allows one to formulate a fusion 
rule of type B for fusion of two points (not associated 
with velocity difference) in the following form: 

F„(ri, r'l; ra, r^; {rfc, r'^}) - F„(ri, r'^; ri, r'j; {rk,r'^^}) 
= 52(|ri -r2|)*„,2(ro;r'i,r'2;{rfc,r',}) cx (;^)^'^^" ■ 

(2.18) 



2. Coalescence of three points (p = 3) 

The next topic of discussion is the fusion of type B 
(without velocity differen ces) of three points, (say ri, r2, 
and ra), as shown in Fig. [IE. As before the separations 
between these points are all of the order of r, and all 
the other separations are much larger, of the order of R, 
R ^ r. As in the case of p = 2 we denote the center of 
mass of these coordinates as rg: rg = (ri + r2 -|- r3)/3. 
We expr ess the two of the velocity differences according 
to (2.16D, and the third velocity difference in analogy: 



(2.19) 



Clearly, the correlation function 



F„(ri, r'l; r2, rs, r'^; {rfc, r'j,}) 

= {w{ri,r[)w{r2,r'2)w{r3,r'3){wr-'') (2.20) 

has three types of contributions. The first, fI^\ is inde- 
pendent of the small separations, 

Ff=F(ro,r;;ro,r^;ro,r^;{rfc,r'J) . (2.21) 

The second type of contribution has either one velocity 
difference or a product of two velocity differences across 
a small separation. Both these contributions have a lead- 
ing term that is proportional to r^^ . The third type of 
contribution has a product of three velocity differences 
across a small separation and is proportional to r^^ . In 
conclusion, the scaling behavior is similar to the one dis- 
played in Eq.(AJ), with the addition of a term constant 
inr, Ff. 



3. General case: fusion of n points without velocity 
difference 

The above discussion of fusion of type B of two and 
three points allows us to offer a general statement con- 
cerning the asymptotic behavior of n-point correlation 
functions Fn when p of the coordinates are separated 
from each other by small distances of the order of r, and 
2n — p coordinates are separated by a large separation of 
the order of R. In terms of the definition (1.6) we will 



take the p close-by coordinates as ri, r2, . . . , rp, and the 
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2n — p coordinates as r-p+i, . . . , r„ and r[, . . . ,r'^. Re- 
call that this choice means that we do not have velocity 
differences across small scales. 

To see the asymptotic beh avior we n eed to repeat 
the substitutions of the type (|2.16| ) and ( 2.19 ), i.e. for 
1<J<P 



'w{rj,rj,) = w{rj,ro) +w{ro,r'^,) . 



(2.22) 



As before the center of mass is Tq = (X]fe=i '''k)/p- The 
result of this substitution can be readily guessed from the 
cases p = 2, 3: 



;r-p,r;;{rfe,rU) = Fl°) 



J=2 



(2.23) 



In this formula i^^^j, is independent of the small sepa- 
ration, and each term F^^^^ originates from j velocity 
differences across small separations. These velocity dif- 
ferences are w(ri,rQ) where i < p. There are always 
Cp = — j)! different products of j velocity differ- 

ences across small scales in F^^p. More explicitly. 



F(2) 

n,p 



A3} ^ 
n,p 



.ii>i2 = l 



^ r. 



.il>i2>i3 = l 



n,3 



(2.24) 
(2.25) 



FIG. [IE The geometry of fusion with a mixture of type A 
and B 

In this Subsection we consider the even more general 
case in which we have within the ball of small radius 
r q pairs of points associated with the velocity differ- 
ences w{ri,r'^ (1 < i < q), and {p — q) points (say 
rg+i, . . . ,rp) associated with velocity d iffere nces across 
large separations of order r, see Fig [IE. The rest of 



the coordinates r/j, r'f., k > p are separated by large dis- 
tances of the order of R, see Fig. 6. For this case we need 
again to reexpress {p — q) velocity differences w{r.j,r'^) 
(for j = q+1, . . . ,p) in the manner of ( ^.22 ), and to sub- 
stitute these into the definition of f „ (l.S). The result 
of this process is an equation similar to ( 2.25 ), but with 
important differences: 

p 

F„(ri,r;;...;rp,r;;{r,.,ra) = E^2p • (2-27) 



There is no constant term now, and the leading order 
contribution is proportional to . The exception (as al- 
ways) is that when q = 1 the leading scaling behavior is 

r^^ . We note that the functions ^ are different from 

the analogous functions i^^^p, but they have the same 

scaling behavior, F^p oc i?*^" (r/ j?)^J . The explicit ex- 
pressions analogous to ( 2.24 2.25| ) can be written down 
when needed. We will find that the subleading terms 
contribute the most important contributions in various 
situations. 



etc. On the LHS of these equations the functions 
depend on all the coordinates. The function Sj is a 
homogeneous function of j separations — ro, with a 
scaling exponent C,j. The functions ^n.j are homoge- 
neous functions of n — p large separations — ro (with 
fc = p + 1, . . . , n and n large separations rj„ — rp (with 
m = 1, . . . , n) with a scaling exponent C„ — Schemat- 
ically we can summarize the scaling behavior of i^^^p as 



F^^' oc 



(2.26) 



We stress that since we a re co ncerned with the limit 
r ^ i?, the leading term in (2.23) is always which is 
independent of r. The leading r dependence is carried by 
F^^'p oc . This does not mean however that the higher 
order terms in (2.23) are unimportant. They will provide 



the leading order contributions to correlation functions of 
velocity gradients as will be shown in subsection 2G. 



E. Mixture of A- and B-type fusions 



F. Fusion rules for the fusion of two or more groups 
of pairs 

The next set of fusion rules is obtained for the struc- 
ture function Fn when two groups of p and q points (with 
p + q < n) tend to and rg respectively. The separation 
between these groups of points is of the order of R. The 
derivation of the fusion rules of type A for the simplest 
situation when all the coordinates are different (and sep- 
arated by r or by R) obviously follows from the same 
basic properties of velocity correlation functions which 
we disc usse d in Subsection 2B. The result looks similar 
toEq. (U): 

F„(ri,r'i;...;r„,r;) (2.28) 
= Fjy{ri, T]^; . . . ; r^, rp)Fq{rp^i, fp^i! . . . ; '^p+g, T^p+q) 

^ ^n,p,q{Xp+q+lit'p+q+l] ■ ■ - ^frnt^n) ■ 

The scaling exponent of ^n,p,q is C n ~ C o — Cq- As in the 



case of the fusion rules (2^), also ( ^.28 ) are not decom- 
positions into products of lower order correlation func- 
tions, and the functions '5' are not correlations of velocity 
differences across large separations. As before the func- 
tions ^ are much larger than the corresponding correla- 
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tion functions in all situations with multi-scaling. Evi- 
dently one can derive similar fusion rules for three, four 
or more groups of coalescing points with large separa- 
tions between the groups. The structure of the resulting 
correlation function will be a product of the correlation 
function associated with each group times some function 
^ of big separations which carries the overall exponent. 

The generalization of these fusion rules (which are of 
type A) to the more complicated cases with fusions of 
type B or to mixed types of fusion is now obvious: we 
can consider every group of point separately in the way 
that we discussed for the case of fusion of just one group 
of points. 



G. Fusion rules for correlation functions including 
unfused velocity gradients 

In this subsection we use the fusion rules obtained 
above to evaluate the leading order scaling behavior 
of correlation functions that include unfused velocity 
derivatives. To be specific, consider the g-order deriva- 
tives ViV2-..Vg, with Vj = d/drj. We are going 
to apply this q-order derivative on correlation functions 
with p fusing points, with q < p, such that the deriva- 
tives operate only on coordinates within this group of p 
points. We will also consider a constrained derivative. 



i.e. such that Vi 



0. With such a 



derivative we get a particularly simple result. There are 
two situations to consider. If the p fusing points undergo 
a fusion of type A, the leading scaling behavior of the q- 
order derivative is simply r^p~i. In the case of fusions of 



type B we consult with Eq.(2.23) and find that the lead- 



ing order contribution with respect to the small distance 
r is 



Vl . . . VqSg{rQ\ri,r2, ■ ■ .,rq)^n,q 



(2.29) 



The contributions arising from the terms F^p all vanish 
under the derivatives. To see that this is so for 2 < j < g 
we recall (cf. ( 2. 24 ), ( 2. 25 )) that although the functions 
i^^-*-* depend on all the p separations — ro, it is a sum 
of functions Sj each of which depends only on a subset of 
j small separations. Thus we find a nonzero contribution 
for the g-order derivative only from terms with j > q. 
The leading contribution for r ^ R always comes from 
the j = q term. It is noteworthy that this contribution 
is independent of the remaining p — q small separations, 
if they exist. 

In summary, the rule is that when we fuse p coordinates 
on which q gradients are applied, the correlation function 
scales asymptotically like r'^p~iJiCn-ip for fusions of type 
A and r'^i-^R'^ri-Qq foj. fusions of type B. The second re- 
sult is independent of the number of additional points in 
the ball of size r that do not have a gradient applied to 
them. These additional points can be even fused together 



or with Tq. This result will be the starting point for the 
discussion of the correlation function jT'n of the balance 
equation. 



III. THE FUSION OF TWO POINTS: TENSOR 
STRUCTURE AND THE EFFECTS OF 
ANISOTROPY 

In the previous section we focused on the scaling ex- 
ponents that characterize the leading contribution to the 
correlation function in the asymptotic regime when some 
points fused together. In this section we address the ten- 
sor structure of the correlation functions, and the sub- 
leading terms that exist because of the dependence on 
the angles between the small separation vector and the 
remaining large separation vectors. This second subject 
is related to the rich and important issue of the decay of 
the effects of anisotropic forcing, and we do not exhaust 
this issue in the present section. Some of the results are 
relevant however in a much broader context. 



A. Tensor structure in the fusion of two points 

In all the previous discussion of the fusion rules when 
two points (say ri and r2) were fused we focused only 
on the scaling exponent of the function F2{p) with 
p = ri — r2. Here we will go further in examining the 
structure of the resulting correlation functions. 

The fusion of two points involves just one small sepa- 
ration distance p = ri — r2- In general the asymptotic 
behavior of the correlation function may depend on the 
angle between p and the remaining large separation vec- 
tors. This dependence is discussed in the next subsection. 
Here we consider the isotropic part which will be shown 
to be the leading contribution. It is easy to determine the 
dependence of F2 on the direction of p using the general 
requirement of incompressibility: 



F^''{p) = A {2 + (2)60.0 -C2 



PaPiS 



nC2 



(3.1) 



This form is standard for the second-order structure func- 
tion in isotropic turbulence. We reiterate that when we 
extract F2 out of a many-point correlation function in 
the process of fusion, there is the issue of the direction 
of p with respect to other vector separations which we 
address next. 



1. Subleading contributions: effect of helicity and anisotropy 

The dependence on the angle of p is an interesting sub- 
ject that deserves full analysis in a separate study. Here 
we only touch on some of the essential issues. 

The analysis of these terms depends very much on 
the nature of the two-point scalar correlation function 
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F2{p) = F^^^p) in anisotropic turbulence. When the 
forcing of turbulence is isotropic, this function depends 
on the magnitude \p\ only. In general however the de- 
pendence on the orientation of p with respect to the 
anisotropic forces may be important. It is useful therefore 
to represent i^2(p) as a "multipole" expansion according 
to 



F2(p) = X]F2,,(p) 



m=-e 



(3.2) 



(3.3) 



In this expansion the p = p/ p, ^ is unit vector and the 
functions l^m are the standard spherical harmonics. We 
chose to expand in these functions since the relevant sym- 
metry group in our problem is the group of rotations 
SO (3). In a scale invariant situation we expect that each 
component F2/ scales like 



F24 oc pf" 



and in general the exponents (3i depend on £. That this is 
so with universal exponents (3i was prove d [|l8[ in the case 
of Kraichnan's model of passive scalar [[19[, but there is 
yet no analogous proof in the case of Navier-Stokes tur- 
bulence. We will assume, in order to proceed, that the 
exponents f3i exist and that they are universal. 

Under these assumption the calculation of the sublead- 
ing contributions to the fusion rules in the case of the 
fusion of two points is straightforward. We first consider 
the partial trace F^" of the n-rank tensor Fn with re- 
spect to the first two indices. Next we decompose it into 
spherical harmonics according to the "multipole" expan- 



Fr(ro + ip,ro-ip;{r,,rU) 

OC 

Y^FZAro + ^p.ro-^p-Ark.r'^}) 



(3.5) 



FZ{ro + kp^ro-^p;{rk,r',})^ ^ 

m=-l 

X /Fr(r-o + |p|,ro-ipeK,r'fe})r,™(Od|, (3.6) 



The first term, -F"^, corresponds to the leading contri- 
bution with the scaling behavior pji?*-''""^^^ which was 
discussed above. (Remember that R is the characteristic 
separation in the correlation function -F""). The next 
order contributions are given by F""^q 



F^",(ri,r;;r2,r^;{r,,r'J) 

= F2,£(p)*„,2,f(ro,;r;,r'2;{rfc,r;,}) 



(3.7) 



The scaling behavior of F2^i can be read from (3.4) im- 
der the usual assumption of universality; the two points 



that fuse together relate to the unfused points in the same 
way that the two-point correlation function relates to the 
anisotropic forcing. We can thus write with impunity: 



F, 



2,1 OC p 



fit 



^ 7: 



,2,e oc 



(3.8) 



To estimate the value of Pi and (32 we need to un- 
derstand what is the physics that determines them. In 
fact, since we considered the partial trace of the corre- 
lation function, we end up with F2 which is even under 
the transformation p —p. Accordingly, although the 
exponents are as stated in (^^) the coefhcients of all odd 
£ quantities are zero. In order to examine odd £ contri- 
butions one needs to form a correlator which is not even 
in p. As an example we consider 



F2"'{P) 



{ua{r + p)up{r) - Ua{r - p)u(3{r)) . (3.9) 



Since this object is manifestly odd in p it vanishes when 
there exists inversion symmetry. For turbulence with 



non-zero helicity F!^^ {p) ^ and the leading contribu- 
(3.4) 

tion to F2^{p) (which is F^f (p)) is determined by the 
flux of helicity. This is reminiscent of the flux of energy 
which determines the leading contribution to the second 
order structure function. Standard K41 dimensional rea- 
soning leads to the prediction /3i = 1, see for example 
|lT| . This holds probably with the same accuracy as the 
K41 prediction for ^2, which is C2 = 2/3 instead of the 
experimental value (,2 ~ 0.70 |l^ , |r^ . 

We can easily determine the tensor structure of 
F^f (p) in isotropic incompressible turbulence (in the ab- 
sence of inversion symmetry) : 



F2.1{p) = ^cp^rPiF^Ap) , F2a{p) oc , (3.10) 

where Cap-y isq the fully antisymmetric unit tensor (ei23 — 
~e2i3 = !)• As we mentioned, in the presence of inver- 
sion symmetr y F 2(—p) = F2(p) and all terms which are 
odd in £ in ( pT^ ) are zero. On the other hand this is 
not the case for F2i{p), which appears in the fusion of 
two points in a many-point correlation function. Even 
if the turbulent flow itself has inversion symmetry, the 
geometry of all the points appearing in (3/7) can lead to 
non- vanishing odd £ components of F2£. The positions 
of the points r'i^r'2 . . . are such that there is no inversion 
symmetry around the center of the fusing coordinates ri 
and r2 which is = (ri -|-r2)/2. Therefore F2e{p) 7^ 
and because of the same constraints it has the same ten- 



sor structure as ( 3.10 ) 



Klip) = ^c.p^p,F2Ap) , F2AP) « p' 



(3.11) 



The physical origin of this term in the multipole expan- 
sion is a local flux of helicity. Even when the average 
helicity flux is zero, the local value of the flux condi- 
tioned on the velocities fixed at certain coordinates may 
be non zero. 



9 



The exponent /?2 is the leading exponent describing the 
rate of decay of the effects of anisotropy, and may be com- 
puted using perturbation theory, disregarding the non- 
perturbative effects which are the subject of this paper, 
see and reference therein. The result is (32 — 4/3, 

and again one expects this resuh to be numerically close 
to the truth. We do not possess presently any numerical 
estimates for (3^ with higher values of i, and as we said 
before it is not guaranteed that these exponents are uni- 
versal. These issues have to be considered independently 
in the future. 



Substituting Eq.( |4.2| ), and considering the stationary 
state in which dFn/dt = we find the balance equa- 
tions 



X>„(ri,r'i; . . .r„,r;j = J'„(ri, r^; . . . r„, r^) , (4.7) 



where the "interaction" term X'n stems from the nonlin- 
ear and pressure terms: 



IV. GENERALIZED BALANCE EQUATION 

A. Derivation of generalized balance equation 

The starting point of this analysis is the Navier-Stokes 
equations for incompressible flows: 



du{r, t) 

dt 



u{r, t) ■ Vu{r, t) + Vp(r, t) = vV'^u{r, t) 



V ■ u{r,t) = 



(4.1) 



In general we need to add a forcing term to these equa- 
tions. It was shown in [Q that as far as the balance 
equations are concerned, the effect of the forcing term is 
felt in the energy containing scales only. For this reason 
we do not write the forcing explicitly. As usual the gra- 
dient of the pressure in (4.1) is eliminated by applying 



the transverse projection operator P. The Navier-Stokes 
equations takes on the form 



du{r, t) 
dt 



P [u{v, t) ■ V] M(r, t) = iyV^u{r, t) . (4.2) 



The application of P to any given vector field a{r) is 
non-local, and has the form: 



Pa(r)], = J dr' P^p{r')ap{r - r'), (4.3) 

where the kernel Pap{r) is the following difference: 

P^r) = 6^pS{r) ^ Pjl^ir) . (4.4) 

Here Papir) is the kernel of the longitudinal projector 
which appears here due to the effect of the pressure term 
in the Navier-Stokes equation: 



.11 



1 



Saf3 



(4.5) 



Given the equation of motion we can take the time deriva- 
tive of Eq.([t|). We find 



dt 



w{ri,r[,t) . 

dw{rj,r'j,t) 
di ■ 



(4.6) 



■'w{rn,r'^,t)) . 



I>„(ri,r;;...r„,r:j =^(-u;(ri,r'i)... (4., 

i=i 



Pu ■ Vu I - l^Pu ■ Vu^ 



■■wirn,r'J) , 



and the "dissipative" term originates from the vis- 
cosity term in the Navier Stokes equation: 



...wir,,r'j)...wirn,r'J) . (4.9) 



In writing these equations we used the fact that in the 
stationary state the time designation is unneeded. We 
also used the following short-hand notation: 



Pu-Vu]^ = / drPa,(3{rj - r)u-,{r)W-^up{r) , (4.10) 



and denot ed b y the Laplacian operator acting on rj 



Equation (4.8) can be written explicitly in the form 



n „ 

P^^"-""(ri,r;;...r„,r:j = ^ / drP^^p{r) (4.11) 
X (u;ai(ri,r'i) . . . Lp(rj,r'j,r) . . . w„„(r„,0} , 



where 



Lf:s{rj,rj,r) = u{rj-r)-'V.jUp{rj-r) (4.12) 



-M(r' — r) • V'jUfj{r'j - r) 



Now we b egin to analyze the balance equation (4.7, 4.£ 
4Jl|, KJ^. 
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B. Galilean invariance of the generalized balance 
equation 

The bal anc e equation must be Galilean invariant. 
Equation (4.9) for J' n depends only on velocity differ- 
ences and clearly is Galilea n inv ariant. This is not so 
obvious in the case of Eq. (4.11) for X>„ because hp in 
( 4.12 ) contains velocities itself. In order to clarify the 
Galilean invariance of T>„ let us express via veloc- 
ity differences only. The first step is to subtract from 
upivj—r) under the derivative Vj the velocity up{r'j —r) 
(independent of rj) and from up^r'^ — r) under the deriva- 
tive Y'j th e velocity ugjr^—r) (independent of r^). Then 
Eq. ([4.1 2D takes the form 



L0irj,rr,r)= [u{r, - r) ■ V , + u{rr - r) ■ V'^] (4.13) 
xwfj{rj - r, r'j - r) . 

For the next step let us introduce u{ri,r[; . . . ; r„, rj^) as 
the mean velocity over all 2n space coordinates: 



k=l 



and the following velocity differences: 



w{rj) = u{rj) 



w{r'A = u{r'A 



(4.14) 



(4.15) 



For brevity we do not display here (and below) argu- 
me nts o f the velocity u. Using ( 4.15| ) we can present 
in ( 4.12| ) as a sum of two terms: 



Lp{r,,rr,r) = L^p'\r,,r'^,r) + L^^\r„rr,r) . (4.16) 
Here the first term depends only on velocity differences: 



[w{rj 



(4.17) 



+w{r'^ — r) ■ V'j]wp{rj — r,r'j — r) . 
However the second term does not have this property: 



(4.18) 



but, as we are going to show, this term gives zero con- 
tribution to Eq . (4.11) for 1>„. Indeed, by substituting 



(4.18) in (4.11) we have: 



n „ 



(4.19) 



X (u;„i(ri,r'i)...L^(rj,r^-,r)...Wc,„(r„,0) . 

In its turn, this equation may be written as the difference 
of two terms, X>^^) = I^^f") - T>^^''\ which correspond to 
the two terms in t he Eq . (4.4) for the kernel Pajpir). 
By substituting in (4.1£) Pa/sir) = 5a(i5[r) we have: 



(2a) _ 



■Wocn {Ttii '''n)) ■ 

(4.20) 

Using the longitudinal projecto r ^^ .^(t') instead of the 



transversal one Pajp{r) in Eq. ( 4.19 ) we have: 

^r^-E/ dr.PH^.ir) (4.21) 

X {wa,{ri,r[) . . . Lf\rj,r'^,r) . . .Wc,„ir„,r'J) . 

Let us show that both these terms are zero (but be- 
cause of different reasons). Consider the first expres- 
sion for 'Dlf"-\ Substituting the explicit form ( [4.18| ) for 
L^^\r j , r'j , 0) and using the incompressibility constraint 
(which allows one to commute Vj + and u) one has: 

n 

^^^=E(^"^(^i'^i) (4.22) 

• ■ • { [i^j + V;-) • u] Wo., (r, , r^.)} ...wo.^ (r„, r'S, 
This equation may be rewritten as 



(4.23) 



Remember that due to space homogeneity the correlation 
function in the second line of this equation is indepen- 
dent of the sum of coordinates Therefore 

'D^^°'^ is indeed equal to zero. 



Consider next Eq. ( |4.2l| ) for X>^^^\ According to 
( 4.18 ) Ljj oc Wf3{rj — r, r'^ — r) and acting on this veloc- 



ity with the longitudinal projector gives zero because of 
the incompressibility constraint. Thus we can conclude 
0. Therefore we can get an ex- 

(1) 



that X», 



■2a) 



-D 



■2b) 



pression for X>„ by replacing L/j in Eq. ( 4.11 ) with L 



taken from Eq. ( [4.171) : 

^ n 

2)r«-"'^(ri,r;;...r„,r;) = - J dr^P„^^(r) 



i=i 



d 

dr. 
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dr 



wp{r.j - r.r'j - r) . . . u;^,. (/•„, r^)^ 



(4.24) 



This expression for X>„ depends only on velocity differ- 
ences and therefore the Galilean invariance becomes ob- 
vious. 
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C. Locality of the Interaction Term 



2. Infrared convergence 



We begin to analyze the interaction terms I?„ (4.24) 
for the most general configuration for which all the 2n 
coordinates r^, are different, and all the n(2n — 1) 
separations are of the same order of magnitude, which 
we designate by i?. This analysis when i? is in the in- 
ertial interval follows the ideas of the analysis presented 
in ^ section 6A for the interaction term of the struc- 
ture function. The interaction term that we face here is 
a significantly more complicated object. The m ain p oint 
is that the integral over r appearing in Eq. (4.24) for 
T^n is "local" in the following sense. First it converges 
in the "ultra-violet" (UV) limit. This limit has to be 
considered when (i) r ^ 0, (ii) when (jj — r) becomes 
very close to any of the 2n — 1 coordinates other than Vj , 
and (iii) when (r^ — r) becomes very close to any of the 
2n — 1 coordinates other than r'y Second, it converges 
in the "infra-red" (IR) limit when r — > oo. The idea for 
the proof of these properties lies in the use of the fusion 
rules which we discussed in Sect. 



1. Ultraviolet convergence 

To demonstrate the convergence of the integral in X>„ 
in the ultraviolet region we can consider any term from 
the sum on j'. Writ ing M^(r, — r) = (X)fe=i "tI'^j ~''"))/'^' 
and using Eq.(4.14), we consider one of the fc-terms in the 
sum. The integral that appears is of the form 



1 



d 



di 



w, 



(4.25) 



w^{rj~r,rk)wp{rj - r,r'j - r) . . . Wa„{rn,r'j'^ . 

As the coordinate r is being integrated over, the most 
dangerous ultraviolet contribution comes from the re- 
gion of small r. In this region the projection operator 
can be evaluated as 1/r^. Other coalescence events of 
r with other coordinates contribute less divergent inte- 
grands since the projection operator does not become 
singular. When r becomes small, there are two possibil- 
ities: (i) rj ^ Tk and (ii) Tj = r^. In the first case the 
correlation function itself is analytic in the region r — > 0, 
and we can expand it in a Taylor series Const + B-r + . . .. 
where B is an r-independent vector. The constant term 
is annihilated by the projection operator. The term lin- 
ear in r vanishes under the dr integration due to r ^ r 
symmetry. The next term which is proportional to is 
convergent in the ultraviolet. In the second case we have 
a velocity difference across t he l ength r. Accordingly we 
need to use the fusion rule (A4), and we learn that the 
leading contribution is proportional to r^^ . This is suf- 
ficient for convergence in the ultraviolet. We note that 
the derivative with respect to rj cannot be evaluated as 
1/r when rj = r^. Rather, it is evaluated as the inverse 
of the distance between rj and the nearest coordinate in 
the correlation function. 



FIG. [VC2 Typical geometry with (n — 1) velocity dif- 
ferences in a ball of radius ii on the left separated by a large 
distance r 3> -R from a pair of points on the right. 

To understand the convergence of X>„ when the inte- 
gration variable r becomes very large we can consider 
again the typi cal term ( [1.25 ). The relevant geometry is 
shown in Fig. IV C 2 . There is one velocity difference 
across the coordinates rj — r and r' — r (which is shown 
on the right of the figure), {n — 1) velocity differences 
across coordinates that are all within a ball of radius 
R (at the left of the figure), and one velocity difference 
across the large distance r which is much larger than R. 
In the notation of this figure the leading order contribu- 
tion for large r is obtained from the fusion rules ( A4 ) for 
the situation on the right and (2.8) for the geometry on 



the left. The resulting evaluation for the leading term is 



r 



C2 



(4.26) 



On the face of it, this term is near dangerous. For K41 
scaling the r-dependence cancels, and the integral is log- 
arithmically divergent. For anomalous scaling the inte- 
gral converges since C»i+i < Cn-i + C2 due to Hoelder 
inequalities. This convergence seems slow. However, the 
situation is in fact much safer. If we take into account 
the precise form of the second-order structure function in 
the fusion rules we find that the divergence with respect 
to rj translates in fact to dSf^'^ {Rj) / dRj^ which is zero 
due to incompressibility. The next order term is conver- 
gent even for simple (K41) scaling. This completes the 
proof of locality of (4.S). The conclusion is that when all 
the separations in £)„ are of the same order of ma gnitude 
i?, the main contribution to the integral in (^^) comes 
from the region r ~ i?. Therefore, the integral can be 
evaluated by straightforward power counting, leading to 



Sn+l{R) 

R 



(4.27) 



It should be stressed that a more detailed analysis 
demonstrates that when the separations p between the 
coordinates that do not involve velocity differences, (i.e 
separations like rjk but not Rj) go to zero, the evalu- 
ation does not change. A direct proof of this fact for 
the case when all such separations are fused (i.e the 
standard structure function) was presented in B. On 
the other hand, if we consider fusions of type A, in 
which coordinates across a velocity difference coalesce 
(p = Rj), we need to be more careful. When two points 
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undergo a fusion of type A the rough evaluation of X>„ 
is [{dS2{p)/dp)][Sn+i{R)/S2{R)] where R is the charac- 
teristic value of large separations. However taking into 
account the tensor structure of the first factor one sees 
that it vanishes due to the incomprcssibility constraint. 
In the next-order term the evaluation of the gradient is 
1/R and X>„ when p ^ 0. 



D. The dissipative term 



The evaluation of the quantity is more straight- 
forward. When all the separations Rj and are of the 



same order i?, the correlator in ( |4.9| ) is evaluated simply 
as Sn{R)- The Laplacian is then of the order of 1/R^. 
The evaluation is 



Sn{R) 
i?2 



(4.28) 



When one of the separations becomes much smaller than 
the rest the evaluation can be read directly from the defi- 
nition (59) and from the fusion rules. Denoting the small- 
est separation by rmin we write 



SnjR) 

^2(i?)rii„ 



(4.29) 



V. THE HOMOGENEOUS EQUATION 

We noted that when all the separations involved in 
our correlation functions are of the same order of mag- 
nitude, and when — > (which is the limit of infi- 
nite Reynolds number Re), the term becomes neg- 
ligible compared to I?„. The ratio JnlT>n is evaluated 
as i^Sn{R) / RSn+i{R), which for fixed R vanishes in the 
limit 1/ — > 0. Thus the "balance equation" becomes a 
homogeneous integro-differential equation 1?„ = which 
may have scale-invariant solutions with anomalous scal- 
ing exponents Cn+i 7^ (n + l)/3. It should be stressed 



that the evaluation (4.27) remains correct for every term 
in 2?„, but various terms cancel to give zero in the homo- 
geneous equation, provided that the scaling exponent C„ 
is chosen correctly. To make this important point clear 
we exemplify it with the simple case n — 2 for which 
X>„ can be greatly simplified. Consider the scalar object 
J^2iri,r[,r2,r'2) = {w{ri\r[) ■ w{r2\r2)). The terms in 
the scalar balance equation for this case are exactly 



V2{ri,r[,r2,r2) 



J2{ri,r[,r2,r'2) 



d[S3iri2') - S3iri2)]/2dri 
d[^3(ri'2) - S3irv2')]/2dr[ , (5.1) 
i^M[S2{ri2') - S2{ri2)] 
V?,[^2(ri'2)-^2(ri'2')]} • (5.2) 



When all the separations are of the order of R we can see 
explicitly that J2 ^ vS2{R)/R^ which is much smaller 



than each term in 'D2. Considering the scale invariant 
solution S^{R) = Ai?^^ where A is a dimensional coeffi- 
cient, we see that 



V2{r 



i,ri,r2,r2) = 



„C3- 

12' 



^C3 — 

' 1'2 



' 1'2' 



Obviously the solution for I?2 = requires the unique 
choice Cs = 1 which is the known exponent for Sz p4[ . 
The coefficient A is now determined as e which is the 
mean energy dissipation per unit mass and unit time. 

We presently do not know how to find the homoge- 
neous solution of the equation = for higher values 
of n. We are even not fully confident that this is an equa- 
tion in the usual sense and not just a constraint that will 
not be sufficient for a unique determination of the scaling 
exponents C„. We feel however that this is an interest- 
ing equation that will offer interesting and worthwhile 
insights. 



VI. THE DISSIPATIVE SCALING FUNCTIONS 

In this section we consider the dissipative "scales" and 
show that they are actually scaling functions. To define 
properly the dissipative length we use the fact that there 
is a cross-over from the scale invariant solution of the 
homogeneous equation to dissipative solutions when J2 
becomes comparable to any of the terms in 2?2- This 
happens when at least one of the separations appear- 
ing in ( ^.2[ ) becomes small enough. Denoting the small- 
est separation as ri„i„ we evaluate J2 ^ 1^5*2 (7'min)/?'min- 
From this we can estimate, using the balance equation, 
S'2(7'min) ~ (S'3 (i?)/i'i?)r^in ^ Ir'^^^Jv. In the inertial 
range we have S2{r) ~ (er)^/'^(r/L)''2~^/'^. The viscous 
scale r\2 for the second-order structure function is then 
determined from finding where these two expressions are 
of the same order of magnitude, i.e. 



lrglv={lr^2fl\rlLf^-^l\ 



(6.1) 



Using the outer velocity scale J7l we estimate e ~ U^/L 
and end up with 



772 ^ LRe 



'1/(2~C2) 



(6.2) 



Note that this result is not in agreement with the ad- 
hoc application of the multifractal model ||lj,|2^,^ which 
predicts 772 - LRe^^/^^+'^^l 

A similar mechanism operates in the general case of 
n ^ 2. As long as all the separations are in the iner- 
tial interval J^n is negligible. When one separation e.g. 
ri2 diminishes towards zero, and all the other separations 
are of the order of R, the internal cancellations leading to 
the homogeneous equation T^n = disappear, and Dn is 
evaluated as in ( |4.27| ). The term is now dominated by 
one contribution that can be written in short-hand nota- 
tion as vV\Fn{ri2, {R})- We can solve for i^„(ri2, {R}) 
in this limit: 
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(6.3) 



(7.7) 



On the other hand we have, from the fusion rule ( ^.28 ), 
the form of the same quantity when ri2 is still in the 
inertial range, i.e. F„(ri2,{i?}) w S2{ri2)Sn{R)/ S2{R)- 
To estimate the viscous scale rjn we find when these two 
evaluations are of the same order. The answer is 



Cn + Cs ~ Cn+1 ~ C2 
2-C2 



(6.4) 



We note that the Hoelder inequalities guarantee that 
Xn > and increases with n. We see that the viscous 
"length" is actually an anomalous scaling function. 



VII. EXACT BRIDGE RELATIONS 

In this section we derive important (and exact) scal- 
ing relations between the exponents C„ of the structure 
functions and exponents involving correlations of the dis- 
sipation field. We consider correlations of the type 



/C(") = {e{x,)w{r,\r[). . MrnK)) oc R 



K 



in) 
It 



{e{xx)e{x2)w{rx\r\). . .i(;(r„|r^)} cx i? ^" 
^ {^t{x,y{x2)...^ix,)w{v,\v\)...w(v,y^)) 



(7.1) 
(7.2) 

(7.3) 



where i? is a typical separation between any pair and 
e(a;) = i^|VM(a;)|^, and we are interested in the scahng 

relations between the exponents \x^^ and the exponents 

(2) 

C„. Note that in this notation is the well studied 
exponent of dissipation fluctuation which is de- 
noted /i. We begin with the rigorous calculation of /ii^^. 

Consider ( p!9| ) for J'„+2(a;i, aj^; 0:2, a;^; ri, r'^ . . .r„,r^J 
in the limit X\ —>■ X2. The leading contribution in the 
limit arises from the Laplacians with respect to the coa- 
lescing points: 

lim Jn+2 = i' hm (V^ + V2)(tt(a;i)M(a;2) 

Xl—^X2 Xi~*X2 

X u;(ri,r;)...t(;(r„,r'J) . (7.4) 

Moving one gradient around and taking the trace with 
respect to the first two tensor indices we see that in this 
limit 



lim JT+2 



-2x;^) 



(7.5) 



As explained in the previous section, when the two points 
coalesce 'Dn+2 of the balance equation loses its internal 
cancellations, and we can therefore conclude immediately 
that 



R 



(7.6) 



In terms of the scaling exponents we are led to the exact 
relation 



(2) 

The scaling relations satisfied by /i„ require consider- 
ations of the second time derivative of the correlation 



(1.8) 



= {w{r^\r[,t)...w{r,\rlt) 
■ ■ ■w{rj\r'^,t) . . .w{rn\r'„,t)). 



(7.1 



Using the Navier-Stokes equations for the time deriva- 
tives we derive a new balance equation X>^^)+B,[f' = Jf ,^^^ 
where, using the definition (4.12), 



/n 
drdr' P{r)P{r')(w{ri\r\) (7.9) 



.L{r„r[,r) . . .L{r.j,r',,r') . . .w{rn\r[. 



Using the fusion rules and following steps similar to those 

described above, we can prove that the integrals over r 

and r' converge. Accordingly, when all the separations 

(2) 

are of the order of i?, every term in Vn is evaluated as 
Sn+2{R)IR^- The term J"^^) takes on the form 



J^^^=-'Y. (V? + V^,)(V| + V,l) (7.10) 

X ('^(^il^'i) • • • i"{r,\r[) . . . w{rj\r'^) . . . u;(r„|0) . 

As before, when all the separation in this quantity are 
of the order of R, the Laplacian operators introduce 
factor of \/R^ and the evaluation of this quantity is 
jT^j^-* ~ v"^ Sn{R) / R"^ ■ Clearly this is negligible compared 

f2l ^9^ 

to typical terms in Vn . The quantity ' contains a 
cross contribution with one Laplacian operator and one 
nonlinear term with a projection operator. The integral 
is again local, and one can show that the evaluation is 
B^^^ ~ vSn+iiR)/ R^ which is also negligible compared 
to typical terms in . 

Now we consider the fusion of two pairs of coordinate, 
e.g. ri2 and — > 0. As before, the cancellations in 
are eliminated, and the evaluation of a typical term 
becomes the evaluation of the quantity. The other two 
terms in the balance equation also become of the same 
order because the Laplacian operators Vf and V3 are 
evaluated as ^-^id '''ai respectively. As before we can 
consider the resulting balance equation as a differential 
equation for i^n (''12 , ^'34 , {R\)- The leading term in this 
equation is 

4z.2v2v2^^„(ri2, r34, {R}) « +I>f ^ Sn+2{R)/R^. 
The solution is 
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^«(ri2,r34,{i?}) rl^rl^Sn+2{R)lv''R^ . (7.11) 

Finally we can write the quantities /C^"-* in terms of the 
correlation function as 

K^-:^ = v^ lim ViV2V3V4:^„+4(ri2,r34,{i?}). 

(7.12) 



Using (7.11) here we end up with the evaluation 



K 



(n) 
'2e 



,n+6 ■ 



(7.13) 



f2) 

For the standard exponent fi = ijq we choose n = 
and obtain the phenomenologically proposed "bridge re- 
lation" 



M - 2 - Ce 



(7.14) 



To our best knowledge this is the first solid derivation of 
this scaling relation. In general, if we have p dissipation 
fields correlated with n velocity differences the scaling 
ex pone nt can be found by considering p time derivatives 
of (1.8), with the final result 



(7.15) 



We see that Eqs.(0, ( p^ and ( [tIsI ) can be guessed if 
we assert that for the sake of scaling purposes the dissipa- 
tion field e(r) can be swapped in the correlation function 
with w^{ri\r[)/ Ri, where i?i is the characteristic scale. 
This reminds one of the Kolmogorov refined similarity 
hypothesis. We should stress that (i) our result does not 
depend on any uncontrolled hypothesis, and (ii) it does 
not imply the correctness of the hypothesis. Our result 
is implied by the refined similarity hypothesis, but not 
vice versa. 

We end this section by noting that the accepted val- 
ues of and ^6 which are about 0.2 and 1.8 respec- 
tively, ar e in g ood agreement with the standard bridge 
relation ( 7.14 ). However, we have presented here a 
"two-dimensional" array of bridge relations depending 
on the indices n and p whose experimental test with high 
Reynolds number flows with good resolution of the dissi- 
pative scales is highly desirable, considering the putative 
exact nature of these relations. 



VIII. SUMMARY AND CONCLUSIONS 

In terms of new predictions the theory described above 
has a lot to offer. Firstly, we have presented the fu- 
sion rules, and it is extremely worthwhile to test them 
against experimental data. One can achieve a reason- 
able test already by using existing data sets from at- 
mospheric boundary layers or grid turbulence. In such 
experiments one measures usually at one space point as 
a function of time. Using the standard Taylor hypoth- 
esis one can measure many-point correlation functions 



for points placed along one line. It is possible to examine 
the properties of such correlation functions when one dis- 
tance is much smaller than all others. Another prediction 
pertains to the viscous scaling function and the anoma- 
lous exponents x„ that characterize them, see Eq.(3.3. 
To test these predictions one needs a good resolution of 
the sub-dissipative scales in a high Reynolds number ex- 
periment. Such data are not readily available, but very 
worthwhile to acquire. Another important point raised 
briefly in this paper has to do with the set of exponents 
(3i which govern the anisotropic properties of the corre- 
lation functions. We noted elsewhere that the same set 
of exponents characterizes the correlation functions 
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APPENDIX A: SOME SPECIAL GEOMETRIES 
OF FUSION AND THEIR IMPLICATIONS 



1. Special geometry of Fig. A 1 



FIG. Al, Special geometry of fusion that is discussed in 



Appendix A 1 



As we discussed in subsection [I C 2 , the evaluation 



(^11 ) is inapplicable when the smallest of the large sep- 
arations (say |ri — r2| = i?min) (which is not associated 
with a velocity difference) becomes similar to small sep- 
arations r (across a velocity difference). To understand 
how to evaluate the correlation function in this case con- 
sider the limit Rmin = w ith the help of the special 
geometry shown in Fig. Al. We have four coordinates, 
ri, r'l, r2 and r2 organized as shown in the figure, i.e. 
r[ — ri along the 2;-axis, r2 = ri, and is on the x-axis, 
with equal distances to ri and r[. This special geome- 
try will help us to derive a result that holds more gen- 
erally. In addition to the velocity differences w{ri,r[) 
and w(r2, r'2) we can have any number of of velocity dif- 
ferences w{ri,r'^), but we demand that the product of 
all the additional n — 2 velocity differences (denoted in 
short-hand as {w}"~^) remains invariant to rotations of 
TT radians around the x axis. In Fig. we show two 
such additional velocity differences across — r'^ and 
r4 — T4. Wc will write the n-th order correlation func- 
tion in short-hand notation as 
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{w{ri,r[)-w{r2,r'^){wr~') 



(Al) 



By definition It) (r2, r'j) = w{r2,r[)+w{r[,r'2). Remem- 
bering that ri — r2 we get 

F„ (r 1 , r i ; r 2 , ; {r fe , r } ) = F„ (r 1 , r'l ; r i , r '2 ; {r fc , rj, } ) 
+F^{ri,r[;ri,r[;{rt,r',}) . (A2) 

Due to the rotation symmet ry the first term on the 
right hand side (RHS) of ( [A^ ) may be written as 
Fnififr'iirijri^; {rk,r'j.}), which by definition is the 
—Fn{ri,r[;r2,r2;{rk,r'i,}). Finally we derive the iden- 
tity 

F„ (r 1 , r'l ; r2 , ; {rfc , r'fe }) = iF„ (r 1 , r'l ; r 1 , r '1 ; {rfc , r'fe }) . 

(A3) 

The object on the RHS has two velocity difi^erences across 
a small separation r, and ?i — 2 large sepa rations. There- 
fore it follows the usual fusion rules ( |2.8| ) for p = 2. Ac- 
cordingly 



Fnirl,r[■,r2,r'2■.,{rk,r'^^}) 
F2{rur[)^nArk,r',}) (x (^)'^' i?^" 



(A4) 



This result was derived for the very specific geometry 
shown in Fig. A 1 . However, the last formula holds much 
more generally. Even if we tilt the vector r[ — ri in an 
arbitrary angle, the result remains invariant. The reason 
is that in the limit r <^ R the anisotropic effects of the 
large scales on the small scales have already disappeared, 
as is shown in Section III. If we ruin the symmetry of ro- 
tation around the x-axis we also do not change the final 
result. The change will be in the factor in Eq.(A3) from 
1/2 to a geometry-dependent factor of the order of unity. 
We can even dissociate r2 from ri over distances of the 



order of \r'i — ri\. Eq.(A4) is rather universal 



2. Special geometry of Fig. 



exists a subleading contribution that is proportional to 
r^3. This is important in the analysis that involves the 
calculation of gradients with respect to these positions. 
To this aim consider the correlation function 

F„(ri,r'i;r2,r^;r3,r^j; {rk,r'^}) (A5) 
= {[w{ri,r'^)-w{r2,r'2)\ [u;(r3, r(j) • i] {u;}"-^) 

where 2; is a unit vector in the z direction. We again use 
a short hand notation {w}""^ for the product of n — 3 ve- 
locity differences across large separations which depend 
on the coordinates r^, r'^ and higher. Using now the fact 
that w{r3,r'^) = w{r3,r[) + w{r'i,r'^) we compute 

F„(ri, r'l; r2, r^; rg, r^; {r^, r'^.}) (A6) 
^{[wir^,r[)-w{r2,r'2)] [w{r[,r'^) ■ £] W"-^) 
+ {[w{ri,r[)-w{r2,r'2)] [w{r3,r[) ■ z]{wr-^) . 

Rotating around the x-axis we can rewrite the last cor- 
relator on the RHS of (A6) as the correlator on the LHS 



with an opposite sign. To see this note that all the terms 
are invariant except the term w{r3,r'i) ■ z that changes 
sign. Thus finally 

Fn{ri,r[;r2, r^; rg, r^; {rfc, r'J) (A7) 
= i ( [wir,,r[) ■ wir2, r'2)] [w{r[,r',) ■ z] {w}-^) 

This correlator has three explicit velocity differences 
across short distances, and therefore according to the 
general rule with j> = 3 it is proportional to r^^ . In 
the general case without rotation symmetry the leading 
term r^^ remains. Therefore we conclude that in the ge- 
ometry in which there are two velocity differences across 
a small separation and one point belonging to a veloc- 
ity difference across a large separation (see Fig. A 2), the 
correlation function Fn can be written to leading order 
as 

i^n(ri,r;;r2,r;;;r3,r^;{rfc,r;,}) = [52(|ri-r2|) (A8) 
+^2(K-r^|) - S2{\ri-r'2\) ^2(K-r2|)] 
*n,2(r'o,r3;{rfc,r;,}) + [S'3(ro|ri, r2, r3) -I- S'3(ro|r'i, r^, r3) 
-S'3(ro|r'i,r2,r3) - S'3(ro|ri, r^, r3)]*„^3(ro, {r^, r^,}) , 



FIG. |A_2[ Special geometry of fusion that is discussed in where = [ri + r[ +r2+r'2 + r^l/b, and ^2, ^3, 



Appendix A 2 



In Fig. 



A 2 we show a situation in which there exist two 



velocity differences across short distances. The difference 
with the general situation shown in Fig. 1 (for p = 2) is 
that we have a fusion of points belonging to short and 
long separations, i. e r2 = r^. The same type of fusion 
existed also in Fig. A 1, but there was only one one short 
distance with a velocity difference across it. Now we 
have two. According to the general rule we should have 
a contribution that is proportional to . The aim of 
this special discussion is to show that in this case there 



and \I'n,3 are homogeneous functions of their arguments 
(in the inertial interval) with scaling exponents C2, Csj 
Cn ~ C2 and — Ca respectively. The function S3 may be 
different from the function S3 in its dependence on the 
angles and the ratios between its argument coordinates. 
But they share the same scaling exponent. 



3. Special geometry of Fig. A 3 
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FIG A i . Special geometry of fusion that is discussed in 



Appendix A 3 

In this subsection we consider the case in which there 
are four coordinates within the ball of size r, but only 
two coordinat es (r i and r[) belong to a velocity differ- 
ence, see Fig. [IE. The other two coordinates (r2 and 
ra) are in the ball, but they relate to velocity differences 
across large separations. To understand the situ ation we 
again consider a special geometry, that of Fig. A3. In 
this geometry ri = r2 = and we study the correlation 
function 

F„(ri,r'i;r2,r^;r3,r?,;{rfc,r;,}) = ([w{ri,r[) ■ z] 
[wir,,r'2)-w{r,,r':,)]{wr~^) . (A9) 

Making the substitutions 

w{r2,r2) = w{ri,r'2) = w{ri,r[) +w{r[,r'2) (AlO) 
w{r3,r'^) = w{ri,r'^) = w{ri,r[) +w{r[,r'3) (All) 

we find that 

2 Fn{ri,r[;r2,r'2;r3,r'^;{rk,r'k}) = {[w{ri,r[) ■ z] 

X |u;(ri,ri)|^{w}""^) + {[w{ri,r\) ■ z] 

xw{r,,r[)-[w{r[,r'2)+w{r[,r',)]{wr-^) . (A12) 

In obtaining this equation we used the symmetry under 
rotation around the x axis in tt. Under this rotation 
ri r[. The first term on the RHS has an explicit 
product of three velocity differences across a small dis- 
tance, and it is therefore proportional to r^^ . The second 
term has two velocity differences across a small scale, and 
according to the previous subsection it contains two con- 
tributions, one proportional to r^^ and the other to r'>-'. 
One can write it in a form similar to (39) and (40), but 
this is not needed at the moment. 
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